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Nucleon drift and diffusion mechanisms in central collisions of asymmetric heavy-ions at near- 
barrier energies are investigated in the framework of a stochastic mean-field approach. Expressions 
for diffusion and drift coefficients for nucleon transfer deduced from the stochastic mean-field ap- 
proach in the semiclassical approximation have similar forms familiar from the phenomenological 
nucleon exchange model. The variance of fragment mass distribution agrees with the empirical for- 
mula cr^^(t) — NcyLc{t). The comparison with the time-dependent Hartree-Fock calculations shows 
that, below barrier energies, the drift coefficient in the semiclassical approximation underestimates 
the mean number of nucleon transfer obtained in the quantal framework. Motion of the window in 
the dinuclear system has a significant effect on the nucleon transfer in asymmetric collisions. 
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I. INTRODUCTION 

In heavy-ion collisions with bombarding energies per 
nucleon in the order of nucleon binding energy, the mean- 
field approach, in terms of time-dependent Hartree-Fock 
equations (TDHF), provides a good approximation for 
describing the average behavior of the collision dynamics 
[ll-Q. For example, the mean- field approximation gives a 
good description of energy dissipation and nucleon drift 
in deep inelastic heavy-ion collisions (DIG). However, the 
collective motion is treated in nearly classical, determin- 
istic manner and the fiuctuations of collective variables 
are severely underestimated [1, Therefore, in the 
mean-field approach, it is not possible to describe en- 
ergy distributions and fragment mass and charge distri- 
butions in DIG. There are other reactions, such as heavy- 
ion fusion at near-barrier energies, spinodal dynamics 
leading to nuclear fragmentation [6|, in which dynam- 
ics of density fluctuations play a dominant role. Much 
work has been done to improve the transport approach 
for describing dynamics of density fluctuations beyond 
the mean-field approximation. Basically, there are two 
different sources of density fiuctuations: (i) fiuctuations 
induces by binary collisions [7- 9] and (ii) mean-field fiuc- 
tuations. Fluctuations and dissipation induced by colli- 
sional mechanism are important at the intermediate and 
high energies. Transport description can be extended in 
a stochastic approach by including the binary collision 
term and its fluctuating part in an analogous manner to 
the Langevin treatment of Brownian motion. In the semi- 
classical limit, this model is known as the Boltzmann- 
Langevin approach. On the other hand, as indicated in 
a recent work [10], the mean-field fluctuations originat- 
ing from the quantal and thermal fluctuations at the ini- 
tial state become the dominant source of density fiuctua- 
tions at low energies. The theory that includes mean-field 
fiuctuations is referred to as the Stochastic Mean-Field 



(SMF) approach. 

In recent works, the SMF approach was employed to 
extract transport coefficients associated with relative mo- 
mentum and nucleon exchange in low-energy heavy-ion 
collisions [ll|,[l3l- Also, some applications of the SMF ap- 
proach have been carried out for analyzing the early de- 
velopment of spinodal instabilities [ij, [ij] . Microscopic 
transport coefficients extracted from the SMF have simi- 
lar forms with those familiar from the phenomenological 
nucleon exchange model [T5l - [l7| , but they provide a more 
refined description of one-body dissipation and the asso- 
ciated fluctuation mechanism. In these initial investiga- 
tions, for simplicity, we calculated transport coefficients 
for symmetric central collisions. It is worth mention- 
ing that the variances of the fragment mass distributions 
calculated with the microscopic diffusion coefficient ex- 
tracted from the SMF approach compares very well with 
the empirical result cr^^(t) = A^exc(i), where a\j^ and 
Ncicc{t) denote the mass variance and the total number 
of nucleon exchange between projectile and target nu- 
clei up to time t, respectively. This observation provides 
another convincing support for the validity of the SMF 
approach. In symmetric collisions, the mass and charge 
numbers of the fragment do not change on the average, 
i.e., the mean value of nucleon drift coefficient vanishes. 

In this work, we investigate nucleon exchange mech- 
anism in central collisions of asymmetric ^°Ga -|- ^°Zr 
system. In Sec. II, we study nucleon drift in the basis of 
the TDHF and compare the results with those obtained 
the semiclassical limit of SMF approach. We illustrate 
the influence of the neck motion on nucleon flux through 
the window area. In Sec. HI, we present a brief descrip- 
tion of the SMF approach and transport coefficients for 
nucleon exchange. In Sec. IV, we investigate the nucleon 
diffusion and variances of fragment mass distributions, 
and conclusions are given in Sec. V. 
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II. NUCLEON DRIFT AND WINDOW 
DYNAMICS 

As an extension of the previous work, we consider cen- 
tral collisions of asymmetric systems at low energies. For 
not too heavy systems, collisions above the Coulomb bar- 
rier lead to fusion, and at below the barrier energies, col- 
liding nuclei exchange a few nucleons and re-separate. In 
order to investigate the dynamics of nucleon exchange, we 
introduce the window between projectile-like and target- 
like nuclei according to the procedure outlined in [lll.ll8|. 
Figure 1 illustrates density profiles at the reaction plane, 
p{x,y,z = 0,t) and the window location in collisions of 
""^Ca -I- ^°Zr system at center-of-mass energy Ecm = 97 
MeV at three different times. The Coulomb barrier is 
Vb = 97.7 MeV. In this figure and in the rest of this 
paper, we perform the numerical calculations with the 
three-dimensional TDHF code developed by P. Bonche 
and his collaborators with the SLy4d Skyrme effective 
force 19]. Since the total nucleon number remains con- 
stant, we can take the mass number of the target-like 
nuclei as the independent variable. The mass number of 
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FIG. 1. Snapshots of the nucleon density profiles on the re- 
action plane, p{x,y,z = 0), are indicated by contour plots 
for the central collision of *°Ca -I- ^'^Zr system at Ecm = 97 
MeV in units of fm~^. The black dot is the center of mass 
point. The red lines indicate the positions of the window xq 
and Do = dxo/dt denotes velocity of the window. 



target-like nuclei is determined by integrating the density 
distribution over the left side of the window, 

AT{t) = J J^^Q{xo-x)f{x,p^,t), (1) 

where xo{t) denotes the position of the neck window. 
The quantity Q{xo — x) denotes the step function and 
f{x,Px,t) is the reduced Wigner function along the col- 
lision direction. Definition of the mass number of target- 
like nucleus is equivalent to 

Arit) = j dxQ{xo- x)p{x,t), (2) 

where p{x,t) = J dp^/ (2Trh)f{x,Px,t) denotes the re- 
duced density of nucleons. 

The Wigner function f{r,p,t) is defined as a partial 
Fourier transform of the single-particle density matrix 
according to 

f{r,p,t)^ J d^sexp (^-lp.s^p(r + |,r-|,i)(3) 

where the single-particle density matrix is defined as 

p{r,r',t) = ^0*,,(r,t)n,(aT)0,,,(r',O. (4) 

iar 

Here, the sum runs over the occupied single-particle wave 
functions with spin-isospin quantum numbers a, r and 
occupation numbers n^. We obtain the reduced Wigner 
function f{x,Px,t) along the collision direction by inte- 
grating over the phase space volume on the window area 
between the colliding nuclei, 

fi^,P.,t) = I dydz^^fir^p,t). (5) 

Using TDHF equations, we can deduce that the rate 
of change of mass number of target-like fragments as 

^^Arit) = VA{t), (6) 

where VA{t) denotes the nucleon drift coefficient deter- 
mined by the net nucleon fiux through the window area, 

f dpx Px - Po f, ,x 

VA{t) ^ - / -^r-^ f(Xo,Px,t). (7) 

J Zirn m 

Here po — mdxo/dt is the velocity of the neck multi- 
plied by nucleon mass. Even though we employ the re- 
duced Wigner function for convenience, the expression of 
the drift coefficient is fully quantum mechanical. Since 
the net nucleon flux is determined by the kinetic term 
in TDHF, Eq. ([7]) does not involve a semiclassical ap- 
proximation. The reduced Wigner function f{x,px,t) is 
exact Wigner transform of the reduced density matrix 
p{x + Sx/2, X — Sx/2, t). The center of mass of the collid- 
ing nuclei in Fig. 1 is located at {Xcm,ycm) = (0,0). The 
red vertical lines indicate the positions of windows xq ■ As 
seen from the figure, the position of the window does not 
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remain constant relative to the center of mass. The direc- 
tions of neck velocities are indicated with arrows during 
the approach and re-separation phase of nuclei. Figure 2 
shows the time dependence of the position of the window 
relative to the center of mass in collisions of '^''Ca + ^"Zr 
system at three different center-of-mass energies. As seen 
from Fig. 2(b), at energies below the Coulomb barrier, 
during the approach phase, the neck moves towards the 
center of mass and then moves away from the center of 
mass. We note that the center-of-mass energy i?cm = 110 
MeV leads to fusion. However, as seen in Fig. 2(a), be- 
fore the system fuses, position of the window undergoes 
a single oscillation. 
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FIG. 2. (color online) The window position is plotted versus 
time for the central collisions of '^°Ca + ^"Zr system at three 
different center-of-mass energies. The center of mass point is 
located at x = 0. 

In order to illustrate the effect of window dynamics on 
nucleon drift, in Fig. 3, we show time evolution of the 
mean mass number of the target-like nuclei. 



AT{t)=AT{0)+ / dsvA{s) 



(8) 



as a function of time at three different center-of-mass 
energies, = 93 MeV (c), = 97 MeV (b) and 

Ecm = 110 MeV (a). In this figure, dotted lines are ob- 
tained by integrating the density over the left side of the 
window from Eq. ([1} and the solid lines are obtained by 
integrating the nucleon flux across the window from Eq. 
([7]). Dashed lines are the result of flux calculations from 
Eq. ([7]), except window motion is neglected. We observe 
that window motion has an important influence on the 
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FIG. 3. (color online) Mean number of nucleon transfer to 
the target nucleus ^°Zr, Arit) — At{0), is plotted versus time 
in central collisions of ^''Ca -|- ^°Zr system at three differ- 
ent center-of-mass energies. The dotted curves are calculated 
directly from Eq. ([T|), the solid lines are obtained by integrat- 
ing the nucleon flux on the right-hand side of Eq. ((7| over 
time, and the dashed lines are calculated similarly with solid 
lines except the motion of the window, second term of the 
right-hand side of Eq. ((TJ, is neglected. 



net nucleon flux or, in terms of transport language, on 
the mass drift coefficient, at both below and above bar- 
rier energies. In this connection, we remark that in the 
phenomenological nucleon exchange model, transport co- 
efficients (mass drift and diffusion) are usually calculated 
by assumin g a fixed window position relative to the cen- 
ter of mass [l6[. As seen in Fig. 3, keeping the window at 
rest relative to the center of mass can introduce an error 
in the order of 50% of the mean number of transferred 
nucleons. 
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A. Nucleon Exchange in the Stochastic Mean-Field 
Theory 

In the SMF approach the mean- field fluctuations 
are incorporated into the dynamical evolution by includ- 
ing the initial state fluctuations in a stochastic approxi- 
mation in a similar manner to the idea presented in Refs. 
[2O-.25] . In the approach, in contrast to the deterministic 
description of the standard mean-field approach, an en- 
semble of single-particle density matrices are generated 
by starting an initial density distribution. A member of 
the ensemble of single-particle density matrices, denoted 
by A, can be expanded in terms of a complete set of 
single-particle wave functions as 



where a and r indicate the spin-isospin quantum num- 
bers of the single-particle wave functions. In the expan- 
sion, elements of density matrices are assumed to be ran- 
dom Gaussian numbers with mean values, 



and variances are determined by 
{6p^^{aT),5p],A^'r')) 



'Srr'Saa' ( 1 " ) -|- (1 



(10) 



Here, nf^ denotes the average single particle occupation 
factors, which are values and 1 at zero temperature 
and given by the Fermi-Dirac distribution at finite tem- 
peratures. Time evolution of the single-particle wave 
functions is determined by the self-consistent mean-field 
h{p^) of the corresponding event. 



d 

ih—(piar{r,t; A) 
ot 



h{p^)(f>,^r{r,t;X). 



(12) 



We introduce the Wigner function f^{r,p,t) corre- 
sponding to each event A of the ensemble of single- 
particle density matrices in the same manner as in 
Eq. ([3]). Also, we define the reduced Wigner function 
f^{x,Px,t) and the mass number of target-like nuclei 
A^{t) in the same manner as in Eqs. ^ and ([T]). The 
time evolution of the mass number of the target-like frag- 
ments in the event A obeys the equation 



dt 



-A' 



dpx Px - Po 
2T:h m 



f^{xQ,Px,t), 



(13) 



which is similar to Eq. ^ except that the nucleon flux 
through the window is calculated by the reduced Wigner 
function in that event. It is possible to convert this 
equation into a Langevin equation for stochastic evolu- 
tion of the mass number of the target-like nuclei. Fluc- 
tuations of the nucleon flux through the window can 
arise from fluctuations of the reduced Wigner function 
originating implicitly from the mass number dependence 



A^^(^) 



^TV-y — ^TV-y ~ Axit) and explicitly from the fluctua- 
tions of single-particle degrees of freedom Sf^{xo,Px,t). 
Then, the reduced Wigner function can be written as 



f^{xo,Px,t) = f{xo,Px,t) + 



f df{xo,px,t) 



OAi 



SA^{t) 



-df\xo,Px,t). 



(14) 



Here, Arit) = (A^it)) and fixo,Px,t) = {f^{xo,Px,t)) 
denote average quantities taken over the generated en- 
semble. For small fluctuations, the ensemble average 
quantities are equivalent to the results obtained by 
ordinary mean-field approximation. As a result, the 
Langevin equation for the nucleon exchange becomes. 



d 



AW) 



VA{t)+l^^yAW)+eAit). (15) 



where va is the drift coefficient for nucleon exchange 
given by Eq. ([7]) in the standard mean-field approach. 
The quantity £,\{t) denotes the fluctuating part of the 
nucleon flux, 



^A 



it) 



dpx Px 



Po 



2TTh 



Sf\xo,Px,t). (16) 



The fluctuating part of the flux is considered as a Marko- 
vian Gaussian random force on the nucleon exchange 
mechanism determined by a correlation function. 



{eAmAit')) = 2Sit-t')DAAit), 



(17) 



where DAA{t) is the diffusion coefficient associated with 
nucleon exchange. In the SMF approach, within the 
semiclassical approximation, it is possible to deduce ex- 
pressions for the diffusion coefficients associated with 
macroscopic variables [ll|, [11] • The expression of dif- 
fusion coefficient for nucleon exchange is given by 



DAAit) = 



dpx \Px - Pol 1 
2Trh m 2 



A+{xo, Px,t), (18) 



where the quantity A^ (xq ,Px,t) is defined as 



A^ixo, Px,t) = ^ I fp'^ixo,Px,t) 
ar ^ 

^^r'i^o-.Px-.-t) 



1 - 



fl{xo,t) 

fp'^{XQ,Px,t) 



ri(xo,i) 



19) 



Here, the mean values of the reduced Wigner distri- 
butions fp^ixQ,Px,t) and f^^ {xQ,Px,i) are defined in 
terms of single-particle wave functions originating from 
projectile and target nuclei, respectively. The quantity 
57(xo,i) denotes the phase-space volume on the window 
and it is calculated with the method explained in Ref. 
fllj . The expression of diffusion coefficient has a sim- 
ilar form with that familiar from the phenomenological 
nucleon exchange model [l6|. However, diffusion coeffi- 
cient extracted from the SMF approach is not restricted 
by adiabatic or diabatic conditions, and provides a more 
refined description of the nucleon diffusion mechanism. 
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FIG. 4. (color online) Nucleon drift coefficients are plotted 
versus time in central collisions of ''"Ca + ^°Zr system at three 
different center-of-mass energies. The solid lines and dashed 
lines are obtained from Eq. ((7| and Eq. pO)) . respectively. 

In the plicnomenological nucleon exchange model, diffu- 
sion coefficient is determined by the sum of fluxes from 
projectile to target and from target to projectile, while 
the nucleon drift coefficient is specified by the net fiux 
through the window. In the spirit of nucleon exchange 
model, we can infer an approximate expression for the 
nucleon drift coefficient as 

^ f^°° \Px-P0\ ,r,n. 

■fAW ~ / TTT ^ {xo,Px,t). (20) 

We note that the expressions for the nucleon drift and 
diffusion coefficients, Eqs. and ([^0]). are valid in 

the semiclassical approximation. In particular, at be- 
low the barrier energies, the reduced Wigner functions 
fp^j,{xQ,Px,t) obtained from the TDHF calculations ex- 
hibit large oscillations in momentum space. In order 
to obtain semiclassical description of transport coeffi- 



FIG. 5. (color online) Mean number of nucleon transfer to 
target nucleus Axit) — At{0) is plotted versus time in central 
collisions of *°Ca -I- ^"Zr system at three different center-of- 
mass energies. The solid lines are obtained by taking time 
integral of the TDHF drift from Eq. ([7|, and the dashed lines 
are calculated from the semiclassical expression Eq. 1)20^ of 
the nucleon drift coefficient. 



cients, we need to smooth out these oscillations in the 
reduced Wigner functions. Similar to the Gaussian over- 
lap approximation introduced in Ref. [l^l, we develop 
a smoothing method of the momentum dependence of 
the reduced Wigner functions, which is presented in the 
Appendix. 



B. Comparison between quantal and semiclassical 
nucleon drift 

Figure 4 shows the drift coefficients as a function of 
time at center-of-mass energies below the barrier, Ecm = 
93 MeV (c) and E^m = 97 MeV (b), and above the 
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barrier energy Ecm = 110 MeV (a). In the figure, the 
quantal drift coefficients calculated from Eq. ([7]) and the 
semiclassical approximation of the drift coefficients cal- 
culated from Eq. (pn| are indicated by solid and dashed 
lines, respectively. The mean-field approach described 
by TDHF equations provides a quantal treatment of the 
single particle motion, and hence we can take the TDHF 
drift coefficient as a reference. As a result, nucleon trans- 
fer via barrier penetration is included into the quantal 
drift coefficients. The semiclassical calculations provide 
a good approximation for the drift coefficients at above 
the barrier energies. However, at below the barrier ener- 
gies, semiclassical calculations underestimate the quantal 
drift coefficients progressively more and more. 

Figure 5 shows comparisons of the mean nucleon trans- 
fer AT{t) — At{Q) calculated from Eq. ([S]) in the standard 
mean-field TDHF and the SMF approaches at the same 
center-of-mass energies below the barrier Ecm = 93 MeV 
(c) and Ecra — 97 MeV (b) , and above the barrier energy 
Ecm = 110 MeV (a). We observe that, at above the bar- 
rier energies, mean value of the nucleon drift is well de- 
scribed by the semiclassical approximation. However, at 
below the barrier energies, quantal effects become impor- 
tant, as a result the semiclassical calculations underesti- 
mate the quantal value of nucleon drift by more than 50% 
a.t E = 97 MeV. At £; = 93 MeV, drift in the semiclas- 
sical calculation goes in the opposite direction. Figure 6 
shows the comparison of asymptotic values of the mean 
numbers of nucleon transfers in the TDHF and SMF cal- 
culations as a function of center-of-mass energy. Below 
barrier energies, the asymptotic values {t — )■ -foo) are 
the mean numbers of nucleon transfer after re-separation. 
On the other hand, at the over barrier energies, for the 
asymptotic values, we take the maximum number of nu- 
cleon transfer at the initial stage of the fusion. At below 
the barrier energies, quantal effects cause the large differ- 
ences in the mean number of nucleon transfers between 
the TDHF and the semiclassical calculations. 



III. DIFFUSION COEFFICIENT AND 
VARIANCE OF FRAGMENT MASS 
DISTRIBUTION 

In this section, we calculate the variance, cr\j^(t) = 
{{A^)^) — (j4y)^, of fragment mass distribution in cen- 
tral collisions of ^°Ca -|- ^°Zr system. It follows from the 
Langevin Eq. p^ . that the variance cr\j^{t) is deter- 
mined by 

J^<A{t) = MtXAit) + 2DAAit), (21) 

where a{t) = dvA{t)/dAT- 

Solid lines in Fig. 7 illustrate the diffusion coefficients 
as a function of time at center-of-mass energies below the 
barrier, E^^ = 93 MeV (c) and E^-^a = 97 MeV (b), and 
above the barrier energy Ecm = 110 MeV (a). The solid 
lines are obtained by Eq. ^TE\\ . while the dashed lines are 
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FIG. 6. (color online) The asymptotic values {t — >■ oo) of the 
mean numbers of nucleon transfer are plotted versus center-of- 
mass energy in collisions of '"'Ca -I- ^"Zr system. The Coulomb 
barrier, Vb, is indicated by the arrow. At energies above 
the barrier, where fusion occurs, the maximum numbers of 
nucleon transfered, for instance the maximum value seen in 
Fig. 5a, are indicated. The solid lines and dashed lines are 
calculated from the standard TDHF Eq. {TJ and the SMF 
approaches Eq. (|20|) . respectively. 

calculated by neglecting motion of the window. We ob- 
serve from this figure that, contrary to behavior of drift 
coefficients, the motion of the window does not affect dif- 
fusion coefficients appreciably. Time dependence of drift 
coefficients VAit) in central collisions of asymmetric ''"Ca 
-I- ^°Zr system at three different center-of-mass energies 
is illustrated in Fig. 4. In particular, at below the bar- 
rier energies, time dependence of drift coefficients in both 
TDHF and SMF approaches arises mainly from time de- 
pendence of the window area and time dependence of 
nucleon density on the window area. Furthermore, be- 
cause of very small value of the mean nucleon transfer, 
mass asymmetry degree of freedom is far from equilib- 
rium. Therefore, the reduced mass dependence of the 
drift coefficients is expected to be very small. In the fol- 
lowing, we neglect the contribution from drift term and 
solve the variance equation ((2T|) to find, 

ai^(t) = 2 / DAAis)ds. (22) 
Jo 

Figure 8 shows variances of fragment mass distribution 
at the same center-of-mass energies below the barrier 
Ecm = 93 MeV (c) and E^m = 97 MeV (b), and above 
the barrier energy Ecm = 110 MeV (a), as a function of 
time. Solid lines in the figure are the result of integration 
over the diffusion coefficients, while solid dots indicate 
the total number of exchange nucleon Ncxdt) until time 
t. The empirical relation cr^^(i) = Ncxc{t) follows from 
the nucleon exchange model and has been widely used 
to analyze the experimental data of deep-inelastic colli- 
sions [2^ [13. Figure 8 illustrates, as in the previous in- 
vestigations for symmetric heavy- ion collisions [l^ . that 
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FIG. 7. (color online) Nucleon diffusion coefficients are plot- 
ted versus time in central collisions of *'^Ca + ^''Zr system 
at three different center-of-mass energies. The solid lines are 
obtained from Eq. (|18|) . while the dashed lines are calculated 
by neglecting motion of the window. 



FIG. 8. (color online) Variances of fragment mass distribu- 
tions are plotted versus time in collisions of ''"Ca + ^"Zr sys- 
tem at three different center-of-mass energies. The dotted 
lines denote total number of exchanged nucleons until a given 
time t, while the solid lines are calculated from Eq. (|22p . 



there is a good agreement between the empirical formulae 
and the calculations based on the SMF approach. Even 
at energies below the Coulomb barrier, the semiclassical 
calculations based on the SMF approach compares well 
with the empirical description iVexc- 



IV. CONCLUSION 

We investigate nucleon drift and diffusion mechanisms 
in central collisions of asymmetric '^^Ca -I- ^"Zr system 
in frameworks of the standard mean-field and the SMF 
approaches. In the SMF approach, we do not carry out 
stochastic mean-field calculations, but we extract diffu- 
sion coefficient associated with nucleon exchange in a 
semiclassical approximation, which has a form similar 



to that familiar from the phenomenological nucleon ex- 
change model. Using this similarity, we also infer an ex- 
pression for the nucleon drift coefficient in the semiclassi- 
cal approximation. In the standard mean-field approach 
provided by the TDHF equations, we can only calculate 
the nucleon drift including quantum mechanical effects 
in the single-particle level, but we can not describe nu- 
cleon diffusion in this framework. At energies below the 
Coulomb barrier, colliding nuclei exchange a few parti- 
cles and re-separate. At above the barrier energies, since 
collisions lead to fusion, we limit our calculations of nu- 
cleon exchange during the early stages of the collisions. 
Our calculations show that not only the size of window 
between nuclei but also motion of window relative to the 
center of mass play an important role in nucleon drift 
mechanism. In the phenomenological nucleon exchange 
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model, the motion of the window relative to the center 
of mass is often neglected. We observe that at below the 
barrier energies, the semiclassical drift coefficient under- 
estimates the nucleon drift deduced from the standard 
mean-field description with TDHF equations. This is 
mainly due to the fact that barrier penetration of the 
single-particles is not included in the semiclassical de- 
scription of nucleon drift. However, in collisions with 
energies above the Coulomb barrier, the semiclassical ex- 
pression for the nucleon drift provide a good approxima- 
tion for the quantal drift deduced from the TDHF frame- 
work. Employing the diffusion coefficient extracted from 
the SMF approach, we calculate the variances of frag- 
ment mass distribution at energies below and above the 
barrier. In such calculations for central collisions, we can- 
not compare the results with experiments, but can com- 
pare them with the empirical relation cr\j^{t) — Ncxc{t)- 
As in the symmetric collisions we investigated earlier 
[Til [T^ , our calculations with semiclassical diffusion co- 
efficients for asymmetric system compare well with the 
empirical relation at energies both below and above the 
Coulomb barrier. This empirical result has been used to 
analyze data in DIC above the barrier energies [l^, [13] . 
However, it is not clear whether such empirical relation 
<^AAi^) — ^cxc(i) is also valid below the Coulomb barrier, 
where quantal effects may become important in nucleon 
exchange mechanism, as we observe in nucleon drift coef- 
ficients. Therefore, in the present calculations, we cannot 
draw any conclusions about the possible quantal effects 
on nucleon diffusion coefficients and the resulting vari- 
ances of fragment mass distributions at energies below 
the Coulomb barrier. 
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Appendix 

We consider the reduced Wigner function 



where 



f^'^{xo,Px,t) = I dydz I ds^exp ( --pxSx 



X 



(l>iar [x - ^,y,Z,tj , (A.l) 

which is associated with the single-particle wave func- 
tions of the target nucleus. It can be rewritten as 

f{xo,p) = j dsexp{-ips/h)F{xo,s), (A. 2) 



z) 



'X(t>i{xQ- -,y,z). (A.3) 

In these expressions for simplicity, we introduce the 
notation f{xo,p) = f^'^{xo,Px,t) and (j)i{xt),y, z) = 
(/)j£rT ixQ,y,z,t) with single-particle wave functions orig- 
inating from the target nucleus. At energies below the 
Coulomb barrier, the reduced Wigner function f{xo,p) 
exhibits rapid oscillations on the window and can take 
negative values due to quantal effects. In order to smooth 
out these oscillations, we need to approximate the quan- 
tity F{xo, s) in terms of a smooth function as a function 
of s, which should be consistent with the following two 
constraints, 

• F*{xq, s) — F{xo, — s), hence the real part must be 
an even function and the imaginary part must be 
odd function. 

• \ims^ooF{xo,s) = 0. 

A possible choice consistent with these requirements is 
to approximate F{xo, s) in terms of a Gaussian function. 
In order to determine the centroid and the width of the 
Gaussian representation, in Eq. (jA.Sp we expand the 
single-particle wave functions around s = up to the 
second order to find, 

Fixo,s)^pixo){l-if3s-^s^). (A.4) 

Here, the quantities a and /3 are given by 

(A.5) 



p{xo) 

1 

Pixo) 



a 



pHxq) +Tx{xo) 



(A.6) 



In these expressions, the reduced mass density, the cur- 
rent density, the squared-momentum density, and the ki- 
netic energy density are defined as 

i 

%{x) = j dydz^^J2{<j)* \7x <j); - <j)i Vx <t>i), (A.8) 

i 

pI{x) dydzY^i^: Vl 0^ + <P^ Vl '/'I),(A.9) 

T^{x)= / dydzJ2\Vx<P^\'^■ (A.IO) 

i 

In mean-field calculations, the expectation value of ki- 
netic energy is computed by using either the reduced 
kinetic energy density 'f^{x) or the reduced squared- 
momentum density p^(a;), since these quantities satisfy a 
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sum rule relation, 



Tx{x)dx. 



(A.ll) 



We note that the reduced kinetic energy density, t^{x), is 
positive for all x values. However, the reduced squared- 
momentum density, p^[x), can take small negative val- 
ues in the vicinity of the window below barrier energies. 
As an example. Fig. 9 shows the quantities t^(x) and 
p^(x) in central collisions of ^''Zr and ^°Ca system at 
a center-of-mass energy of Ecm — 97 MeV and _at time 
t = 10.5 X 10"^^ s. We observe that the quantity p^^j;) for 
target as well as projectile can take small negative value 
in the vicinity of window position xq. Consequently, the 
parameter a can take negative values. 



^ 1.5 



0) 

T3 
>> 
CD 

_o 

g 
1^ 



0.5 




10 15 



FIG. 9. Quantities (solid line) and Tx{x) (dashed line) 

are plotted versus position x in collisions of ''"Ca + ^"Zr sys- 
tem at _Bcm = 97 MeV and at time t = 10.5 x 10~^^ s for 
target and projectile. The neck position xq is indicated by an 
arrow on the figure. 

In the Gaussian approximation, in order to produce 
the correct second order Taylor expansion, we introduce 
a new parameter 7 = 0; — 2/3^. When a > and 7 > 0, 
we can approximate the quantity _F'(xo,s) in terms of a 
single Gaussian as 



F{xa,s) « p{xo) exp {-il3s - ^s^) , 



(A.12) 



By taking the Fourier transform, we obtain a Gaussian 
form for the reduced Wigner function as well. 



f{xo,p) ^ 2p(xo)^/-exp 



(p-m 



21 



(A.13) 



When the quantity a < 0, we can write the expansion 
Eq. (|A.4|) as 



F{xo,s) w p{xo) [l-i/3s+ —s' 



(A.14) 



There are many possibilities to approximate F{xo,s) in 
terms of Gaussian functions that have the same second 



order Taylor expansion with Eq. (|A.14I) and that are 
consistent with the two constraints indicated above. We 
consider a double Gaussian approximation according to 



F(xo, s) = p{xo) exp {-il3s) 



l7l 



2exp(-^s2)-exp(-^s2) 



which leads to the Wigner function given by 



(,A.15) 



f{xo,p) = pixo) 




4 exp 



■ exp 



(A.16) 



Figure 10 shows an illustration of the reduced Wigner 
functions given by Eqs. (jA.13[) and (jA.16[) by employing 
the same /3 and the same magnitude of 7. We observe 
that for 7 < 0, the reduced Wigner function exhibits 
a negative tail, which is closely related to the quantal 
behavior of the Wigner function. In the calculations pre- 
sented in this work, relatively small negative tails in the 
reduced Wigner functions do not make a significant ef- 
fect on the semiclassical expressions of the nucleon dif- 
fusion and drift coefficients given by Eqs. and (pOj) . 
Substituting Eq. (|A.13I) and/or (|A.16|) for the reduced 



m 

CO 
i— 



10 
8 
6 
4 
2 

-2 



. ' 1 ' 


1 1 1 1 1 1 1 1 1 1 1 1 

/A Y<0 - 




/ \ 7>0 ■ 




1/ ^\ 
A ^\ 

/ / \ ^ — 
/ \ 


, 1 . , , 1 , , . 1 , . , 1 . 


-1 


1 2 




p (arbitrary units) 



FIG. 10. Reduced Wigner functions, Eq. (|A.13|) and Eq. 
()A.16|) . are plotted in arbitrary units with parameters 'p{xo) = 
1, /3 = 0.5, and 7 = ±0.5. 

Wigner functions originating from target and projectile 
into Eq. ([T8|. we obtain the nucleon diffusion coefficient 
as 

Daa =T{pp,[3p,'yp) +T{prp,/3T,jT) 

+n(7p,7T), (A.17) 

for 7p > and jt > 0, 

i^^^ = 2r(pp,/3p, |7p|/2) - r(pp,/?p,2|7p|) 
-|-F(pT,,/3T,7T) + 2n(|7p|/2,7T) 
-n(2|7p|,7T), (A.18) 
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for 7p < and 7^ > 0, and 

I?AA = 2r(pp,/3p,|7p|/2)-r(pp,/3p,2|7p|) 
+2r(p^,/3T, |7t|/2) - r(pr,/3T,2|7T|) 
+4n(|7p|/2,|7T|/2) + n(2|7p|,2|7T|) 
-2n(|7p|/2, 2|7t|) - 2n(2|7p|, |7t|/2), (A.19) 

for 7p < and 7t < 0. The functions T and 11 are 
defined as 



r(p,/3,7) 



n(7p,7y) 



Zm 



exp 



7 



/3erf 



V7 



(A.20) 



o ^ '^^P 

TO S2 7P + 7T V 7^^ 7t 



v^(7-p/3t + 7T/3p) 
1 H , exp 



X erf 



7p/3t 



'7t) 
■ 7T/3p 



(7p/3t + 7t/3/ 



\ \/"/pjt{ip + it) ^ 
where the error function is given by 



erf(a;) = I e 



it) 
(A.21) 



'dy. 



Employing the formula given in Ref.jll| for the phase 
space volume VL(xQ^t) on the window, we obtain for fol- 
lowing expressions for the nucleon drift coefficient, Eq. 

(uni), 



VA = T{pp,(3p,-/p) - V{pj,,I3t-,1t) (A.23) 

for 7p > and 7^ > 0, 

VA = 2r(pp, /3p, |7p|/2) - r(pp, /3p, 2|7P|) 

-r(py,/3T,7T) (A.24) 

for 7p < and 77- > 0, and 

VA = 2r(pp,/?p, |7p|/2) - r(pp,/?p,2|7p|) 

-2r(pj,,/3T,|7T|/2)+r(pj,,/3T,2|7T|)(A.25) 



(A.22) for 7p < and 7^ < 0. 
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